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Simple model for QCD analysis of the proton helicity structure 
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In this paper we use the experimental data to obtain the polarized parton distribution functions 
(PPDFs) in the LO and NLO approximations. The analysis is based on the Jacobi polynomials ex- 
pansion of the polarized structure function (PSF). Our calculations for polarized parton distribution 
functions based on the Jacobi polynomials method are in good agreement with the other theoretical 
models. 
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I. INTRODUCTION 

The theoretical and experimental status on the 
spin structure of the nucleon has been discussed 
in great detail in several recent reviews (see, e.g., 
Refs. P, H, 0] ■ Deeply inelastic scattering provides 
a clean way to extract the parton densities of nucle- 
ons. After the initial observation that the nucleon 
spin is not formed by the quarks dominantly Q, de- 
tailed measurements of the polarized structure func- 
tions followed during the last 20 years. The central 
question concerns now the parton distribution func- 
tions and their scale evolution rather than just their 
first moment. Since the nucleon spin receives also con- 
tributions from the angular momentum of the quarks 
and gluons, these degrees of freedom have also to be 
studied. During the recent years several comprehen- 
sive analysis of the polarized deep inelastic scatter- 
ing (DIS) data, based on next-to-leading-order quan- 
tum chromodynamics have performed. In these anal- 
ysis the polarized parton density functions are either 
written in terms of the well-known parameterizations 
of the unpolarized PDFs or parameterized indepen- 
dently, and the unknown parameters are determined 
by fitting the polarized DIS data. 

Determination of parton distributions in a nucleon 
in the framework of quantum chromodynamics (QCD) 
always involves some mo del- dependent procedure. In- 
stead of relying on mathematical simplicity as a guide, 
we take a viewpoint in which the physical picture of 
the nucleon structure is emphasized. That is, we con- 
sider the model for the nucleon which is compatible 
with the description of the bound state problem in 
terms of three constituent quarks. We adopt the view 
that these constituent quarks in the scattering prob- 
lems should be regarded as the valence quark clusters 
rather than point-like objects. They have been re- 
ferred to as valons. The idea of nucleon as a bound 
state of three quarks was presented for the first time 
in Ref . Q . In the valon model, the proton consists of 
two "up" and one "down" valons. These valons thus, 
carry the quantum numbers of the respective valence 
quarks. Hwa and et al. Q found evidence for the 
valons in the deep inelastic neutrino scattering data, 
suggested their existence and applied it to a variety 
of phenomena. In fs*] unpolarized PDFs and hadronic 



structure functions in the NLO approximation were 
extracted. In Ref. @ the polarized valon model is 
applied to determine the quark helicity distributions 
and polarized proton structure functions in the NLO 
approximation by using the Bernstein polynomial ap- 
proach. The extraction of the quark helicity distri- 
butions is one of the main tasks of the semi-inclusive 
deep inelastic scatterin g (S IDIS) experiments (HER- 
MES COMPASS [nl, SMC [121) with the po- 
larized beam and target. Recently in Ref. [H, |3| 
the polarized valon model was applied and analyzed 
the flavor-broken light sea quark helicity distributions 
with the help of a Pauli-blocking ansatz. The re- 
ported results of this paper are in good agreement 
with the HERMES experimental data for the quark 
helicity distributions in the nucleon for up, down, and 
strange quarks from semi-inclusive deep-inelastic scat- 
tering [lO|. 

Since very recently experimental data are available 
from the HERMES collaboration ;i5S| for the spin 
structure function gi, therefore there is enough mo- 
tivation to study and utilize the spin structure and 
quark helicity distributions extracted via the phe- 
nomcnological model. We performed the first recal- 
culation for the spin structure function gi u p to NLO 
[Toj using the Jacobi polynomials expansion 17 1. This 
paper is based on the Ref. [l^ . 

The plan of the paper is to give a brief review of 
the theoretical background of the QCD analysis in 
two loops. The method of the QCD analysis of polar- 
ized structure function, based on Jacobi polynomials 
are written down in Section 3. A description of the 
procedure of the QCD fit of gi data and results are 
illustrated in Section 4. 



II. THE THEORETICAL BACKGROUND 

Let us define the Mellin moments for the polarized 
structure function g\{x, Q^)'. 



x''-'glix,Q^)dx 



(1) 



The contribution to the structure function gi{N, Q^) 
up to NLO can be represented in terms of the po- 
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larized parton densities and the coefficient functions 
ACf" in the Mellin -N space by [1] 



form for gi(N,Q^) involves some unknown param- 
eters. If the parameters can be obtained then the 
computation of all moments of the PPDFs and the 
polarized structure function, gi{N,Q'^), are possible. 



Ag(A^,g')] + ^2AC'; Ag{N,Q')} , 

(2) 

in this equation the NLO running coupling constant 
is given by 



1 



/3iln(lnQVA|^) 

/3olnQVA|j^ /33(lnQVA|^)2 



(3) 



The symbol Ag denotes the strong coupling constant 
normalized to Ag — as/(47r). Notice that in the above 
the numerical expressions for /3o, Pi are 



Po = 11 - 0.6667/, 
Pi = 102 - 12.6667/ , 



(4) 



where / denotes the number of active flavors. In 
our calculation , we choose Qo = I GeV^ as a 
fixed parameter and A is an unknown parameter 
which can be obtained by fitting to experimental data. 

In Eq. 0, Aq{N, Q^) = Aq^N, Q^) + Aq{N, Q^), 
Aq{N, Q^) and Ag{N, Q^) are moments of the polar- 
ized parton distributions in a proton. AC^ , ACg 
are also the iV-th moments of spin-dependent Wilson 
coefficients given in Ref . [l^ . 

According to improved polarized valon model 
framework, determination of the moments of parton 
distributions in a proton can be done strictly through 
the moments of the polarized valon distributions. 

The moments of PPDFs are denoted respec- 
tively by: Au„(iV,g2), Ad^{N,Q'^), AE(A^,Q2) and 
Ag{N,Q^). Therefore, the moments of the polarized 
u, d, S and g density functions in a proton are: 

Au^iN, Q2) ^ 2AM^^(iV, g2) X AA%/p(Ar), (5) 



Ad,{N, Q2) = AM^^(7V, Q2) X AM^/^iN), (6) 



AE(iV,Q2) = AM^iN,Q'){2AAC^j,iN)+AM';,/p{N)), 



(7) 



Ag{N,Q') = AMS«(Ar,Q2)(2AMy/p(iV)+AAf^/^(iV)) 



(8) 

In the above equations M'^j^{N) and Mj^^{N) are 
the moments of polarized valon distributions, which 
introduced in Ref. [191 . It is obvious that the final 



III. THE METHOD OF THE QCD ANALYSIS 

One of the simplest and fastest possibilities in the 
PSF reconstruction from the QCD predictions for its 
Mcllin moments is Jacobi polynomials expansion. The 
Jacobi polynomials are especially suited for this pur- 
pose since they allow one to factor out an essential 
part of the x-dependence of the structure function into 
the weight function [3]. Thus, given the Jacobi mo- 
ments a„((5^), a structure function f(x, Q^) may be 
reconstructed in a form of the series [19l|-[23l| 

xfix, Q2) = _ ^ a„(g2)e-.'5(x), (9) 

71=0 

where N^ax is the number of polynomials and 9"'''(x) 
are the Jacobi polynomials of order n, 



(10) 



where cj"'' (a, P) are the coefficients that expressed 
through F-functions and satisfy the orthogonality re- 
lation with the weight x^{l — x)" as following 



dx x^ii ~ xret'^{x)ef-'^{x) = Skj , (11) 



For the moments, we note that the dependence is 
entirely contained in the Jacobi moments 

a„(g') = / dxxf{x,Q^)Q'^-^{x) 



= Y.cf\a,P)f{3+2,Q^), (12) 

3=0 

obtained by inverting Eq.( [9]), using Eqs. (flOl 
[TT|) and also definition of moments, /(j, Q^) = 
/g dx x^~^xf{x, Q^). 

Using Eqs. ([9l fT2|) now, one can relate the PSF with 
its Mellin moments 

Xfff— (X,g2) = X'^(l-X)" J2 ^n^(x)x 

71=0 



Y,cf\a,P)gi{j + 2,Q^), 



(13) 
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FIG. 1: NLO polarized parton distributions at the input scale Qq^I.O GeV^ compared to results obtained by BB model 
(dashed line) (ISET=3)[3l, AAC (dashed-dotted line) (ISET=3)i39j, GRSV (dashed-dotted dotted line) (ISET=1)[43| 
and LSS (dashed-dashed dotted line) (ISET=l)l4H. 



where f/i(j+2, Q^) are the moments of polarized struc- 
ture function. N^axj ol and /3 have to be chosen so 
as to achieve the fastest convergence of the series on 
the R.H.S. of Eq. (fT3l) and to reconstruct xg\ with the 
required accuracy. In our analysis we use N^ax = 9j 
a = 3.0 and /3 = 0.5. The same method has been 
applied to calculate the nonsinglet structure function 
xF^ from their moments [24-31]. 

Obviously the Q^-dependence of the polarized 
structure function is defined by the Q^-dependence 
of the moments. 



IV. THE PROCEDURE OF THE QCD FIT 
AND RESULTS 

For the QCD analysis presented in Ref. the fol- 
lowing data sets were used: the HERMES proton data 
[H, [3pL the SMC proton data jllK the E143 proton 
data j34| , the EMC proton data j35l. [36| . In the fit- 
ting procedure, using the CERN subroutine MINUIT 
[37| . we defined a global for all the experimental 



data points and found an acceptable fit with minimum 
xVd-O.f . = 0.978 in the LO case and xVd.o.f . = 0.933 
in the NLO case. In Table. 2 of Ref. lj| we presented 
and compared the results which are based on Bern- 
stein and Jacobi approaches. 

In Figures 1 the parton distribution functions in 
next-to-lea ding order for all sets of parameterizations 
[38, 3l M 

l4lj and their errors are presented at the 
starting scale Qq. 
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